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Abstract: The methods for 
obtaining the gain pattern of an 
Air Traffic Control Radar using 
the data collected by NUSAT 1 
are presented. The basic 
equations for this are given; 
these require an independently 
determined satellite attitude. 
In addition we show how 
simultaneous measurements made 
by separate channels can be 
related to obtain a single axis 
attitude, i.e., the direction of 
propagation of the incident 
radiation, and the power density 
of the radiation. This is done 
by solving a system of three 
algebraic equations. 

Introduction . NUSAT 1 was the 
first satellite to be launched 
from a Getaway Special 
Cannister. It has six L-band 
channels which are driven by six 
orthogonally directed antennas. 
The purpose of these systems is 
to detect Air Traffic Control 
Radar (ATC) transmissions and 
measure the power density of the 


incident radiation, thus allowing 
the gain pattern of the radar 
transmitter to be determined. 
Since the satellite is free 
floating (no attitude control) 
and since the gains of the 
satellite's receiving antennas 
are not unidirectional a 
fundamental problem arises in the 
determination of the direction of 
propagation (with respect to a 
satellite based coordinate 
system) of the incident 
radiation. To surmount this 
problem a crude attitude 
determination system was included 
in NUSAT 1. It will not give an 
attitude accurate enough for data 
reduction . 

The purpose of this paper is 
twofold; one, we will present the 
equations needed to calculate the 
gain of the radar transmitter and 
power density of the radar wave. 
These can be used if an 
independently determined attitude 
is available. Two, we will also 
show how three simultaneous 
measurements made by three 
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separate channels can be related 
to obtain a single axis 
attitude; i.e., the direction of 
propagation of the incident 
radiation, and the power density 
of that radiation. This is done 
by solving a system of three 
algebraic equations. In 

addition, we will give a short 
error discussion to indicate the 
type of accuracy available. 

1 . Overview of NUSAT 1 For 

those readers who are unaware of 
the need and current methods for 
obtaining the power pattern of 
an ATC radar, we will give a 
brief summary of the situation; 
[1] contains more detailed 
information . 

Air Traffic Control Radar 
The ATC radar transmitting 
antennas rotate at a rate of one 
revolution in four to ten 
seconds. They have an elevation 
beam width of about 40 degrees 
(the lower 15-20 degrees is used 
to detect aircraft) and an 
azimuth beam width of 2 degrees. 
They transmit on a frequency of 
1030 Mhz . If the antenna is 
pointing too low, then ground 
interference causes nulls and 
lobes. So for safe operation in 
controlling air traffic, the 
radiation pattern for many 
different azimuthal orientations 
of the radar antenna is desired. 
In lieu of this, at least the 
location of the nulls and lobes. 

The Satellite . The shape of 
NUSAT 1 is a twenty-six sided 
polyhedron. Upon launch, the 
satellite was not given any 
initial spin, nor is there any 
attitude control. Thus, the 

satellite's attitude is allowed 
to drift . 

The satellite has a rather 
crude attitude determination 
system - eight symmetrically 
located wide angle field of view 
sensors (detecting radiation 
primarily in the visible region 
of the spectrum). 



Flo. 1. Partial po war gains in tha horizontal Q* and 

vtrtlcal Q« . 


The features of present 
interest to us are the six L- 
band antennas whose major axes 
are directed in six orthogonal 
directions. Each is a circularly 
polarized double Archimedes 
spiral antenna with a nominal 
total gain of 1.4 dB . One cut of 
the pattern (partial power gains 
in the horizontal and vertical 
directions) of a typical antenna 
is depicted in Fig. 1. Note that 
the usable range is a hemisphere 
centered about the antenna's 
major axis and that the gain is 
down by 6 or 7 dB at 60 degrees. 
2. The Basic Equations. The 

calculation used to determine the 
gain of a radar antenna can be 
split into two steps. First, the 
measured value of the received 
power is used to find the power 
density of the incident wave at 
the satellite. Then this value 
is combined with data from 
orbital prediction routines to 
compute the gain of the radar. 
These calculations are based on 
the transmission equation. 

The fundamental equation 

governing power transmission is 
the Friss Transmission Equation 
[2, Sec 2.17] : 

P r =(A/4 J rr) 2 G r (^ r ,^ r )x 

G t (tf t .<MP eP t (1) 

where 

P r - power supplied by the 
receiving antenna, 
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P t *= power applied to the 
transmitting antenna, 
r = distance between the 
antennas , 

e — total efficiency, 
p - polarization loss 
factor , 

G r (0 r ,<£ r ) — isotropic gain of 

the receiving antenna in 
the (8 r ,<j> v ) direction, 

G t^t>^t) “ isotropic gain of 
the transmitting antenna in 
the direction. 

The relation between these 
parameters is shown in Fig. 2. 

The power density W of the 
incident electromagnetic wave is 
related to the received power by 

W=P r 4*-A- 2 (G r (0 r ,^ r )p) | - 1 . (2) 

By using this in equation (1) , 
one can show that the desired 
quantity G t (0 t ,^ t ), the gain of 
the radar antenna - ground 

interference in the direction 
(8 , is 

G t(0t>^t> = ^Trr 2 W/P t e (3) 

In logarithmic decibel form this 
becomes 

[G t ( 6 t , <j> t ) ] dB = 201ogr + [W] dB 

+ 10 log ( 4 tt / e P t ) . 

where the last term may be 
dropped if the absolute gain is 
not desired and the quanities 
Ft> e are nominally constant. 

In our application, P r is 

measured by the satellite; then 
equation (2) is used to 
calculate W. In (2) the angles 
8 x ,<f> r represent the direction of 
propagation of incident wave; 
i.e., the direction of the ATC 

radar in a satellite based 

coordinate system. These can be 
found from an attitude 

determination system or they can 
be computed simultaneously with 
W. In addition, the 

polarization loss factor, p, 



Antenna 

Fig. 2. Relation of parameters 
for power transmission. 

must be evaluated. The 

evaluation of p and computation 
of W will be discussed in more 
detail in the next section. Once 
W is found, it can be used in (3) 
to get the gain of the radar 
antenna. Here, the variables 

r > 0 t > ^ t are the spherical 

coordinates of the satellite in a 
local coordinate system based at 
the radar antenna. These can be 
found from orbital prediction 
routines . 

3 . ... Computation of power density 
arid the polarization loss factor. 
The power density W of the 
incident radar wave is computed 
using (2). Nominally, the radar 
wave is linearly polarized and 
the receiving antennas are 
circularly polarized; so p may be 
approximated by 1/2. For more 
accuracy p needs to be evaluated; 
this may be done using the gain 
measurements performed on the 
receiving antennas. 

The antennas were measured 
using a modified version of the 
Mu 1 t ip le -Ampl itude - Component 
method. The entire gain pattern 
(letting 9 x ,<f> r range over useful 
values, a hemisphere about the 
major axis) of each antenna was 
mreasured with a linearly 
polarized transmitter oriented at 
angles i/> of 0, 45, 90, and 135 

degrees (see Fig. 3). This 

resulted in a measurement of the 
partial power gains , G G q , 
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G .^25 respectively in the 

direction (0 r , <£ r ) . 



propagating from direction 


( ^ r » ^ r ) • 

Letting L, in Fig. 3, 

represent the electric field of 
the incident wave allows the 
polarization ratio to be related 
to either of the angles r/> or 5. 
By using the angle and, 

assuming that the wave is 

perfectly linear; i.e., that the 
power in the direction 

perpendicular to L is zero, one 
can show that 

9 9 9 

p — cos ^ +r z, s in^ + 

rsin2^cosj0/(l + r^) 

where r , /? are the polarization 
ratio and relative phase angle 
of the receiving antenna in the 
direction (0 r> ^ r ). A similar 
equation is available for p 
using the angle 8. Note that 
for a perfectly circular antenna 
we have cos/? = 0 and r= 1 so that p 
reduces to 1/2. We also note 
that there is the option of 
appr oximateing cos/? with a 
nominal value for each antenna. 
If this is done, only the 
measurements , Gq need to be 
performed and 

G r p = G^cos^i/) + Ggsin^V’ 

+ r s in2V>co s/3G r / ( 1 + r ^ ) . 


Since a complete set of four 
measurements are available for 
NUSAT 1, both r and can be 

found for each direction (0 r> ^ r ). 
Omitting the intermediate steps, 
we get two different equations to 
use for the product G r (0 r ,<£ r )p in 
equation (2). Using the angle ^ 
we get 

G-rP = G,icos 2^ + 6 0 sin^+ 

<f/2)si„2j(G 45 -G 135 ). < 4 > 

Using the angle <5 we get 

G rP = G 45 cos ^ +G i 

(5) 

+(l/2)sin2S(G0 -G^) . 

If the satellite's attitude can 
be found, then 0 r ,</> r can be 
determined from the relative 
location of the radar and the 
satellite at the time the 
mesasurement was made. Further, 
V> and 8 can be found by applying 
the attitude matrix to a vector 
in the direction of the local 
vertical at the radar antenna. 
Thus, G r(^r»^r) can be determined 
from either equation (4) or (5). 
For better accuracy, the equation 
that has the smallest value for 
sin2^ or sin26 should be used 
(subtraction produces a loss of 
significance in the third terms). 

Alternately, the power density, 
W, and the angles 0 r ,</» r can be 
determined simultaneously by 
relating three received power 
levels. For this though, the 

polarization loss factors must be 
approximated by 1/2. Recall that 
NUSAT 1 has six L-band antennas 
oriented so that their major axes 
are directed in the six 
orthogonal directions. When the 
satellite is not in a null of the 
radar's radiation pattern at 
least three channels will be able 
to detect a signal. The geometry 
is shown in Fig. 4, where the 
major axes of the antennas which 
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Fig. 4. Three L-band antennas and the direction of the 
incident wave with power density w. 

are able to detect a particular 


s ignal 

have been placed a 

long 

the x , 

y , z directions . 

From 

equation ( 2 ) we have for 
three channels : 

the 

p i - 

(A 2 /4tf)Gi(0 ,^) Pl W 


p 2 = 

(A 2 /4jt)G 2 (0 ,^)p 2 W 


p 3 = 

(A 2 /4tt)G 3 (0 ,^)p 3 W 


where 

® i > ^ i » P i » i — 1,2,3 

are 


associated with channel i and 
9 ,4> give the direction of the 
radar in body coordinates. If 
the polarization loss factors 
are approximated by 1/2 this 
system can be solved for W , 6 , <f> . 
The gain functions G^(9,<f>) are 
only known through test 

measurements, yet any of several 
numerical methods can be used to 
solve this system. If the 

chosen method requires the 
partial derivatives of the 
functions Gi , then it might be 
necessary to smooth these 
measurements. Actually this 

system should be solved in 
logarithmic decibel form: 

[ p i]dB = 101og(A 2 /4ir) 

+ [ Gj_ ( 0 , ^) ] dB + [ p i ] d g 

+ [ w ] dB > i = 1 .2,3. 

Combining each t?i]dB with 
1 Olo g ( A 2 /4»r ) and using to 

denote the combined quantity 
this system becomes, with the dB 
subscripts dropped 


Pi = Gi(9 ,<f>) + W + p 1( 

j> 2 _ G 2 (6,<f>) + W + p 2 , ( 6 ) 

P 3 = G 3 (0,<£) + W + p 3 . 

For actual computation, each p 
must be approximated by -3dB. 

4. Preliminary error discussion . 
We now give an error discussion 
to indicate the origin of major 
errors and to realize the type of 
accuracy available. We will only 
give a short discussion of the 
errors as the complete error 
analysis is quite lengthy and is 
still in the process of being 
completed . 

The gain of the radar antenna 
in the direction (6 t ,<f> t ) is 
computed from (3) with the power 
density being calculated from 
(2). From (3), the relative 

error in G t (9 is: 

A G t /G t = 2Ar/r + AW/W + Ae/e 
+ A p t/ p f 

The angles and the 

distance r are found from orbital 
prediction routines which are 
usually quite accurate. The 

efficiency e and power supplied 
to the transmitter P t are close 
enough to being constant that 
they will also not contribute 
much error. The relative error 
in the power density AW/W will 
be at least two orders of 
magnitude larger than the others. 
It can be found from (2), 

AW/W = AP r /P r + 2AA/A + 

AG r p/G r p . 

the received power P r is obtained 
from measurements of the signal 
in the L-band channels . A 

calibration can be carried out 
for each channel, thus keeping 
AP r /P r comparatively small. The 
relative error in the wave length 
AA/A will be very small, the 
relative error in neglecting the 
doppler shift is .01 %. The term 
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G r (0 r ,<£ r ) will produce the 

largest error. 

An error analysis for 
equations (4) or (5) is in 

progress, but no results are 

available at this time. An 

indication of the error was 


obtained by considering G r and p 
separately. If the angles 0 r ,<f> r 
are found independently from W, 
then the error in G r (0 r ,<£ r ) 

will depend on the errors in 0 r 
and 4> r ; the accuracy of the 
gain measurements of the 

receiving antennas; and the 
variation of G with respect to 
0 r ,<£ r . If the angles 0 r , <£ r are 
known to be within .5 degrees, 
then AG r /G r will be bounded by 
12% since the gain measurements 
indicate that the variation is 
bounded by .25 dB/deg. Some 
quick estimates of the error in 
approximating p by 1/2 indicate 
that it will average about 1.2dB 
but may be as high as 2 dB . 

Thus, the total error in 
finding the radar gain G t (0 t ,<£ t ) 
given values of 0 r ,<£ r accurate 
to within .5 degrees and 
aproximating p by .5 will be on 
the order of 1 . 6 dB but could be 
as much as 2.3 dB . 

To give a preliminary 
indication of the error in 
determining 0 r ,<£ r and W 

simultaneously, the author 

worked several examples. These 
gave the following bounds: 

AS , A(f> < 10 degrees, 

AW < ldB. 

5 . Conclusion As yet the 

satellite has not made any 
actual radar measurements. 

Therefore, no test results are 
available for presentation. 
But, as reproduced here, the 
necessary equations are in place 
for calculating the ATC radar 
gain in case some measurements 
are made. These same equations 
will be needed for the second 
satellite to attempt an ATC 
radar calibration, NUSAT2. 


1. C.A. Bonsall and R.G. 
Moore, NUSAT 1 - The first 
ejectable Getaway Special, 
Proceedings of the Getaway 
Special Exper imenter ; s Symposium, 
Goddard Space Flight Center, 
October , 1985. 

2. C. Balanis, "Antenna Theory 
Analysis and Design," Harper and 
Row , 1982 

3. "IEEE Stanard Test Procedures 
for Antennas: IEEE Std 149- 
1979, Published by IEEE, Inc. 
1979, Distributed by Wiley 
Inter s c ience . 


198 


